Electronic transport through a Ceo molecular bridge: The role of single and multiple 

contacts 
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The effects of different contact geometries, bond dimerization, and gate voltage on quantum 
transport through a Ceo molecule are studied by the Landauer-Biittiker formula based on the Green's 
function technique. It is shown that the number of contact points between the device electrodes and 
the molecule can play an important role in the electron conduction. The transmission is due to the 
resonant tunneling when the electrodes are contacted to one carbon atom of the molecule. In the 
case of multiple contacts, the interference effects are responsible for the change of the transmission 
through the Ceo- The bond dimerization and a gate voltage shift the molecular levels and by 
adjusting the related parameters the electron conduction can be controlled. 



I. INTRODUCTION 

Molecular electronics using single molecules as active 
elements is a promising technological concept with fast 
growing intereslji'^. Recent improvements in manipula- 
tion of individual or small numbers of molecules permit 
us to contact molecules to metallic electrodes and mea- 
sure their electronic transport properties^'^'^'i. Elec- 
tronic transport through single molecules strongly de- 
pends on the nature and quality of the contacts with 
electrodes. If, e.g., a molecule is weakly coupled to the 
electrodes, the charge at the molecule becomes strongly 
localized and the transport takes place in the regime of 
Coulomb blockade. In the opposite case (strong cou- 
pling), we can expect to approach the ballistic regime. 

The selection of a molecular bridge and the accurate 
control on its coupling to the electrodes are basic prereq- 
uisites for designing and manufacturing single molecule 
electronic devices. In such a structure, the energies of 
the molecular orbitals, in particular, the highest occu- 
pied molecular orbital and the lowest unoccupied molec- 
ular orbital (LUMO), are of crucial importance for the 
electronic transport through single organic molecules. 
Among many types of molecules, the fuUerene Cgo is suit- 
able for molecular bridge because its LUMO is situated 
at relatively lower energies in comparison with the other 
organic molecules^. The electronic structure of isolated 
Ceo molecules shows a large gap (« 2 eV), which indi- 
cates that the molecules should behave like insulators at 
room temperature^. However, when they are contacted 
by metallic electrodes, charge transfer occurs and they 
become conducting through the LUMO of the isolated 
molecules. 

During the last decade, the electron conduc- 
tion through a Cgo molecule has been extensively 
investigated both experimentall y^"'^^ '^'^'^'^'^ 
and theoretically^ i^'''^^i^^i^°i^^'^^ . Joachim and co- 
workersiSiii studied the conductance through a Ceo 
molecule sandwiched between the Au(llO) surface 
and a scanning tunneling microscopy (STM) tip at 



ambient temperature. The results showed that the 
current-voltage [I-V] characteristics of the molecule are 
linear at low voltages due to the absence of molecular or- 
bitals around Fermi energy. Porath and co-worker o^^i^^ 
deposited isolated Ceo molecules onto a gold substrate, 
covered by a thin insulating layer. In this way, a double 
barrier tunnel junction configuration was realized in 
which a Ceo molecule is coupled via two junctions to 
the gold substrate and the tip of a STM. The tunneling 
/- V spectra exhibited a nonvanishing gap in the curves 
around zero bias due to the single-electron tunneling 
effects, such as the Coulomb blockade and the Coulomb 
staircase. Park et al.^'^ fabricated single- molecule tran- 
sistors based on individual Ceo molecules connected to 
gold electrodes. The device showed strongly suppressed 
conductance near zero bias voltage followed by steplike 
current jumps at higher voltages. Such a gap could 
be reduced to zero by adjusting gate voltage. More 
recently, Neel and co-workers^^'^® studied contacts to 
a Ceo molecule on Cu(lOO) and observed that the 
conductance rapidly increases in the transition region 
from tunneling to contact, with decreasing tip- molecule 
distance. 

On the other hand, using the tight binding model, 
it was theoretically found that the strength of the 
metal/Ceo interaction and the geometry of the contact 
between the tip and the molecule play an important role 
in the drastic increase in the conductance^. Based on 
the Green's function method and the Landauer-Biittiker 
formula, it has been shown that a loop current emerges 
in a Ceo molecule when the electron energy approaches 
the energy levels of the molecule. The magnitude of the 
such loop currents can be much larger than that of the 
source-drain current^. By incorporating an extra atom 
at the center of the fuUerene molecule, it is possible to 
control the currents in the loops and hence the procedure 
of transpori^i^. By using density functional theory, the 
electronic transport through a Ceo molecule in between 
carbon nanotube leads^^ and Al metallic electrodes has 
also been investigated22i2ii22. 

The Ceo molecule consists of 12 pentagons and 20 
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FIG. 1: The three difFerent ways of coupling between the 
Ceo molecule and the one-dimensional electrodes used in this 
work. Black circles show the position and the number of cou- 
plings. 

hexagons. Due to dimerization, the carbon-carbon bonds 
in the molecule have different lengths: ri = 1.46 A for 
the single bonds (bonds on the pentagons) and r2 = 1.40 
A for the double bonds (bonds on the hexagons that are 
not shared by a pentagon). We believe that the effect 
of bond dimerization on the electron transmission may 
considerably affect the I-V characteristics, under suit- 
able conditions. Furthermore, when two electrodes are 
connected to the molecule, the number of contact points 
will be dependent on the direction of the molecule; thus, 
single or multiple contacts may occur. Therefore, the ef- 
fects of contact geometry in the presence of bond dimer- 
ization and a gate voltage on the electron conduction 
through the molecule should be studied. Such features 
have not been investigated in the above-mentioned stud- 
ies. It is the purpose of this paper to study the role of 
bond dimerization, multiple contacts (see Fig. 1), and 
gate voltage in the coherent quantum transport through 
a Cgo molecule, based on the tight-binding model and 
the nonequilibrium Green's function technique. 

The paper is organized as follows. The theoretical 
model and formalism are given in section II. In sec- 
tion III, we present the numerical results of the coherent 
transport through the Cgo molecular bridge with differ- 
ent electrode/molecule contacts in the presence of bias 
and gate voltages. A brief conclusion is given in section 
IV. 



II. MODEL AND METHOD 

We consider a system consists of a Ceo molecule 
attached to one-dimensional metallic electrodes. The 
Hamiltonian of the system is described within the tight- 
binding approximation with only one orbital per atom, 

H = eic\c, - Ujcjcj . (1) 

i <ij> 

Here, is the on-site energy and will be set to zero ex- 
cept in the gated region (Cgo molecule) where it is equal 



to Vg- cJct (cicr) is the creation (annihilation) operator 
for an electron at zth site and tij is the hopping matrix 
element between nearest- neighboring sites i and j. The 
hopping strength in Cgo molecule depends on the C-C 
bond length; thus, we assume different hopping matrix 
elements: ti for the single bonds and t2 for the double 
bonds. In the case of bond dimerization, it is reasonable 
to use t2 — (ri/r2)^ ti^^. The coupling between the near- 
est sites in the electrodes is taken to be t = ti and that 
between the molecule and the electrodes is taken to be t'. 
In this study, we assume that the electrons freely prop- 
agate and the only resistance arising from the contacts. 
This means that the transport is ballistic^; therefore, we 
set t' = t/2 according to Ref4. 

The Green's function of the Cgo molecule coupled to 
the two metallic electrodes (source and drain) in the pres- 
ence of the bias voltage is given as 

Gc(e,K) - [el-Hc-±Lie-eVa/2)^±nie+eVa/2)]-' , 

(2) 

where He describes the Hamiltonian of the molecule in 
the absence of the electrodes and T,l and E/j describe the 
self-energy matrices which contain the information of the 
electronic structure of the electrodes and their coupling 
to the molecule. These can be expressed as 

'^L,r{() ^ TCL,RgL,R{f)TL,RC , (3) 

where gL,R are the surface Green's functions of the un- 
coupled leads, i.e., the left and right semi- infinite leads, 
f is a matrix that couples the molecule to the leads and is 
determined by the geometry of the molecule-lead bond. 
Note that, in the semi- infinite one-dimensional electrodes 
described by the single-band tight-binding model, only 
the first site is connected to the molecule. As a result, 
the surface Green's function for the semi-infinite leads 
can be written as gL,R = — (1/t) e''^^'^''', where A:i(_R) is 
the wave vector in the left (right) electrode^i. 

When the molecule is brought close to an electrode, the 
bonding between them will depend on the molecule ori- 
entation. This orientation can be such that only one car- 
bon atom, a pentagon or a hexagon, of the Ceo molecule 
be in contact with the leads. Therefore, one can expect 
different conduction through the molecule, which arises 
due to the interference effects, and will be discussed in 
detail in the next section. For this reason, we make use of 
the nonequilibrium Green's function technique to obtain 
the current for a constant bias voltage Va between two 
electrodes. Since the total Hamiltonian [Eq. ([T])] does 
not contain inelastic scatterings, the current is computed 
from the Landauer formul a^^i^^ 

2p 1"°° 

I=-f^J Tie,Va)[f{e-eVa/2)-f{e + eVa/2)]de, (4) 

where / is the Fermi function and T(e, Va) is the energy- 
and voltage-dependent transmission function given by 

T(e,K) = Tr[fL(e-eK/2)Gc(e,K)ffl(e+eK/2)Gj,(e,K)] • 

(5) 
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FIG. 2: Transmission probability as a function of energy at 
Vg = 0.0 V and Va = 0.5 V for different contacts in two cases: 
t2 = 1.1 ti (solid curve) and t2 — ti (dashed curve). 
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FIG. 3: The same as Fig. 2 but for Va = 2.0 V. 



III. RESULTS AND DISCUSSION 



The coupling matrices T^^n, also known as the broad- 
ening functions, are related to the self-energies through 
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(6) 



Equations ^ and ^ form the basis for our analysis 
of the coherent transport through the molecular bridge. 
Our approach, as a real-space method, makes it possible 
to model arbitrary the number of contacts. In this re- 
gard, the core of the problem lies in the calculation of 
the self-energies f^L,R- In the case of contact through 
a single carbon atom of the molecule, only one element 
of the self-energy matrices is non-zero. However, for the 
transport through opposite pentagons or hexagons, 25 
or 36 elements of the self-energy matrices are nonzero, 
respectively. 



The present formalism can be applied to the systems 
in which an arbitrary number of carbon atoms of the Ceo 
molecule can be coupled to the one-dimensional metallic 
electrodes. In this study, the coupling through one, five, 
and six carbon atoms to the electrodes will be considered 
(see Fig. 1). These cases were chosen as the most proba- 
ble experimental orientations^-^. In the numerical calcu- 
lations we set ti = 2.5eVi2, t2 = 1.1 ti^^, Ep = 0.0 eV, 
andT = 300K. 

Connecting the molecule to the electrodes broadens 
and shifts the discrete states of the molecule. This broad- 
ening of the states depends on the strength of the cou- 
pling to the contacts as well as the wave function of the 
particular state. Therefore, one can expect different be- 
haviors for the transmission coefficient (TC) of the sys- 
tem in the cases of single and multiple contacts. In order 
to investigate such a behavior, in Figs. 2 and 3 we have 
shown the TC as a function of the energy of the elec- 
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tron which is emitted from the left electrode into the 
molecule, for K = 0.5 V at Vq = 0.0 V and Vq = 2.0 V, 
respectively. It is clear that the transmission functions 
have large values (peaks) near the molecular levels of 
CgQ. In the coherent transport, the electron wave func- 
tion is assumed to extend coherently across the whole 
system. Accordingly, when the electron energy nearly co- 
incides with the molecular levels, the electron resonantly 
transmits through the molecule and a large transmission 
arises. From Fig. 3, it is evident that, a gate voltage 
shifts the states and, due to this effect, the transmission 
channels can significantly vary . Therefore, a gate po- 
tential can change the electron conduction through this 
device, producing field-induce molecular switching. 

The number of peaks in the case of connection to five 
carbon atoms of the Cgo molecule is lower with respect to 
the other cases. In the case of connection to one carbon 
atom, the value of TC at the Fermi energy is near zero, 
however, it is considerable in the case of five contacts. In 
both figures, the main factor for the difference between 
the single and multiple contacts is the interference ef- 
fect. In fact, when the molecule is contacted through one 
carbon atom to the electrodes, the transmission through 
the molecule corresponds to the resonant tunneling ef- 
fects. With increasing the number of contact points, the 
interference effects around these contact points become 
important, some resonances might completely disappear, 
and the transmission curve changes. By comparing the 
Figs. 2(b) and 2(c) with Fig. 2(a), and Fig. 3(b) and 
3(c) with Fig. 3(a), one can easily observe the influence of 
interference effects on the electronic transport. The phys- 
ical meaning of the interference effect is that the electron 
waves in the molecule which come from the different con- 
tact points may suffer a phase shift. Thus, a constructive 
or destructive interference in the propagation process of 
the electron through the molecule can occur. As men- 
tioned in Sec. II, the effect of contacts is described by 
the self-energy matrices. Therefore, the Green's function 
(and hence the density of states) of the coupled molecule 
and the transmission spectrum vary with the number of 
contact points^"*. 

Note that in the quantum transport, the width 
of the transmission resonances depends on the elec- 
trode/molecule coupling strength. In one-dimensional 
electrodes, only one atom (end site) is coupled to the 
carbon atoms. Therefore, if we increase the number of 
atoms at the surface of electrodes, by choosing electrodes 
with finite cross sections, then the width of resonances 
becomes significantly broader for the cases of five and six 
atomic contacts^*'. We have also studied the effects of 
changing the electrode/molecule coupling strength. The 
results for the case of single contact showed that the TC 
increases with increasing t' . For strong coupling, the 
peaks are broadened which indicate that the electronic 
transport can no longer be considered as resonance tun- 
neling through eigenstates of the isolated Cgo molecule^. 
In the case of multiple contacts, however, the broaden- 
ing of peaks and the increase in TC are not noticeable. 
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FIG. 4: Current vs voltage characteristics at Vg ~ 0.0 V for 
different contacts in two cases: t2 = 1.1 ti (circle) and t2 = ti 
(square) . 



This suggests that the interference effects play a domi- 
nant role in the transmission through the molecule, when 
we increase the number of contact points. 

On the other hand, it is important to note that the 
presence of both single and double bonds causes a shift 
in the peaks of TC. Shift to the higher or lower energies 
depends on the peak positions. From the Figs. 2(a) 
and 2(c), it is clear that the difference in bond lengths 
gives rise to an extra peak in the TC curves due to the 
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FIG. 5: The same as Fig. 4 but for Vg = 2.0 V. 



appearance of a sharp dip at the energy near -2.63 eV. 
The gate voltage shifts the position of peak, as shown 
in Figs. 3(a) and 3(c). Therefore, the effect of bond 
dimerization (i.e., <2=l-l^i) is that, the degeneracy of 
one of the molecular levels is lifted^. By applying a 
suitable gate voltage, this effect may be important, if 
this level be nearly coincided with the Fermi energy of 
the metallic electrodes. 

In order to provide a deep understanding of the elec- 
tronic transport, we have shown the I-V characteristics 
in Figs. 4 and 5, at Vg = 0.0 V and Vg = 2.0 V, respec- 
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FIG. 6: Gurrent as a function of gate voltage at different 
contacts in the case of t2 — 1.1 fi. 



tively. Since our structure is symmetric, we obtained a 
symmetrical behavior in the I-V curves with respect to 
Va = 0.0 V. For this reason, we have not shown the results 
for negative applied voltages. At Vg = 0.0 V and in the 
case of one contact point, a threshold voltage is needed to 
generate current through the device which shows that, at 
low applied voltage the device is in its off state. For a cer- 
tain gate voltage such as Vg = 2.0 V, the device is turned 
on and the current linearly increases for low applied volt- 
ages [see Fig. 5(a)]. For the case that t2 = 1.1 ii, the 
threshold voltage is nearly two times larger than that in 
the case of t2 = ti. When the number of carbon atoms 
which take part in coupling between the molecule and 
the electrodes increases, the I-V characteristics show an 
Ohmic behavior at low applied voltages. Such a behavior 
is reasonable because in this situation, the hybridization 
with the electrodes is stronger and there are more paths 
for the electrons to pass from the metal to the molecule^^. 
The curves in all figures, particularly in Fig. 4, show a 
steplike behavior which indicates that a new channel is 
opened. 

The present results for the I-V characteristics of the 
cases of single and multiple contacts are qualitatively in 
agreement with the experimental results in Refs.— i^^, 
andlS respectively. In Ref. ^, the bond length difference 
has been ignored because they believed that the shift in 
the energy spectrum is less than 2%, and there is no 
significant difference in the electron conduction through 
the molecule. The present results, however, show that 
by adjusting the parameters, a considerable difference in 
the I-V characteristics between the cases of <2=^i and 
t2=l.lti may occur. 

In Fig. 6, we have shown the current as a function of 
gate voltage. The results reveal that the electron conduc- 
tion can be strongly dependent on the gate voltage. For 
the values of the gate voltages that the current becomes 
zero, the molecule behaves as a semiconductor and for the 
other values it acts as a metal. Such a behavior can be 
clearly seen in the cases that one or five carbon atoms of 
the molecule are contacted to the electrodes. Therefore, 
the results suggest that the Ceo molecule is an interesting 
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candidate for operation of devices as a nanoscale current 
switch. 



IV. CONCLUSIONS 

Using the nonequihbrium Green's function technique 
and the Landauer theory, we have investigated the effects 
of different contact geometries on the electronic trans- 
port through a Ceo molecule sandwiched between two 
metallic electrodes. It has been shown that the trans- 



mission curves and the I-V characteristics calculated at 
the single and multiple couplings to the electrodes can be 
completely different due to the interference effects. In ad- 
dition, the influence of gate voltage on Ceo and the effect 
of different bond lengths in this molecule were observed 
as a shift in the transmission peaks. Our results therefore 
indicate that the contact geometries, bond dimcrization, 
and charge-transfer doping with gate voltage play impor- 
tant roles and may change the physical picture of electron 
conduction in Ceo-based molecular devices. 
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